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Abstract 

We consider a general scheme for calculating the leptonic constant of heavy quarkonium 
QQ in the framework of nonrelativistic quantum chromo dynamics, NRQCD, operating as the 
effective theory of nonrelativistic heavy quarks. We explore the approach of static potential 
in QCD, which takes into account both the evolution of effective charge in the three-loop 
approximation and the linearly raising potential term, which provides the quark confinement. 
The leptonic constants of bb and cc systems are evaluated by making use of two-loop anoma- 
lous dimension for the current of nonrelativistic quarks, where the factor for the normalization 
of matrix element is introduced in order to preserve the renormalization group invariance of 
estimates. 
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1 Introduction 

A description of leptonic constants for a heavy quarkonium composed of heavy quark and heavy 
anti-quark QQ plays an important role in studies of QCD dynamics in the sector of heavy quarks. 
So, the precision calculations of heavy quark masses m^ c and running coupling constant a s in 
QCD are performed in the framework of QCD sum rules H] and effective theory of nonrelativistic 
heavy quarks, NRQCD j2|. In such the approach the experimental data on both the masses of 
S'-wave quarkonium states and their leptonic constants |§ are used. The consideration of NRQCD 
sum rules was done for the bottomonium with taking into account of one-loop gluon corrections 
in [|J (next-to-leading-order, NLO) as well as with two-loop corrections in [0-0] (next-to-next-to- 
leading-order, N 2 LO), and for the charmonium in |S|]. In general, the sum rules of NRQCD could 
be used for the solution of inverse problem, i.e. the calculation of leptonic constants at given 
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values of heavy quark masses, a definite normalization of coupling constant and the experimental 
mass spectrum of quarkonia. The only additional entry is caused by that the analysis in the 
framework of NRQCD sum rules supposes that the region under study remains correctly defined 
only if one takes into account the higher excitations of ground state, so that one should fix the 
ratios of leptonic constants, i.e. their relative weights in the sum rules, along with the mass 
spectrum. Another way is to go in the region of parameters, where the contributions of ground 
state excitations are negligible, and the term of gluon condensate becomes significant, while its 
Wilson coefficient requires cumbersome calculations in higher orders of perturbation theory. 

An alternative method is the description of heavy quarkonium in the framework of potential 
approach. In this way, if the static potential in QCD with infinitely heavy sources in the funda- 
mental representation of color gauge group is given, then in the leading order of NRQCD, i.e., 
if we neglect relativistic corrections and interactions depending on the quark spins, the heavy 
quark masses can be determined with a high accuracy from the comparison of experimental mass 
spectrum of heavy quarkonia bb and cc with the estimates in the NRQCD. In this procedure, one 
calculates the spectra of stationary energy levels E for the nonrelativistic Schrodinger equation 
with the static potential, so that the masses are defined by the expression M = mi + m 2 + E, 
where are the masses of heavy quark and heavy anti-quark composing the quarkonium. In 
addition, in ref. || authors derived the static potential in QCD, so that it takes into account recent 
perturbative results at short distances 110, |HJ as well as the regime of confinement in the infrared 
region. This potential is straightforward modification of model by Buchmiiller and Tye fl2| , who 
suggested the method with unified two-loop /^-function for the effective charge a. This function 
satisfied two definite asymptotic limits at a — > (the perturbation theory) and a — > oo (the po- 
tential confining the quarks in hadrons, in the form of linearly raising term under the distance 
increase). The calculations in three loops [[[(], [Tl| shown that the value and sign of coefficient 02 
modelled in the paper by Buchmiiller and Tye are significantly different from the correct results for 
the /3-function. The modification allows us to reach a consistent description for combined data on 
the mass spectra of heavy quarkonia and the evolution parameter Aqcd for the running coupling 
constant measured at high virtuality of rriz ||. Thus, in the potential approach we have got a 
possibility to estimate the leptonic constants of heavy quarkonia in NRQCD with the quark masses 
and coupling constant a s adjusted with a high accuracy. 

In this way, we consider a general scheme for the calculation of leptonic constants, since the 
results of NRQCD have to be related with the quantities defined in full QCD. There are two 
challenges in this respect. The first one is the Wilson coefficient depending on the scale of calcu- 
lations in NRQCD, which relates the current of nonrelativistic quarks considered in the leading 
order of inverse heavy quark mass, with the electromagnetic current of heavy quarks in full QCD. 
An anomalous dimension of this coefficient /C is known with the two- loop accuracy [ X3— Jl5j . The 



factor /C is matched with the value, which is given by the equality of currents in QCD and its 
nonrelativistic approximation at a scale /ihard about the mass of heavy quark. Therefore, the 
operator equality of currents in full QCD and NRQCD is completely defined up to the two-loop 
accuracy. The second challenge is the matrix element for the current of nonrelativistic quarks 
over the vacuum and physical state of quarkonium. Because of renormalization group invariance 
of physical quantities, this matrix element has got the anomalous dimension equal to that of the 
current in NRQCD. For the first time, in the potential approach we perform the consideration 
involving the factor A providing the renormalization invariance of results. In this way, we have to 
introduce a scale of normalization for A, so that the reasonable region of scale change is restricted 
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by the physical characteristics for the quark system under study. Then, the obtained estimates of 
leptonic constants for the quarkonia are stable in a region of change for both the QCD-to-NRQCD 
matching scale /ihard and the point of perturbative calculations in NRQCD /if ac t- 

In Section 2 we consider the description of leptonic constants for the heavy quarkonia in the 
framework of NRQCD. The potential approach is actually given in Section 3. The numerical 
estimates are presented in Section 4. In Conclusion we summarize the obtained results. 



2 Leptonic constants 

In the NRQCD approximation for the heavy quarks, the calculation of leptonic constant for the 
heavy quarkonium with the two-loop accuracy requires the matching of NRQCD currents with the 
currents in full QCD, 

J? CD = Q7.Q, J, NRQCD = 

where we have introduced the following notations: Q is the relativistic quark field, x an d 4> are the 
nonrelativistic spinors of anti-quark and quark, = o v — v u (a ■ v), where v is the four- velocity of 
heavy quarkonium, so that 

J? CD = /C(/i hard ; /i fact ) • J u NRQCD (^ ct ), (1) 

where the scale /ihard gives the normalization point for the matching of NRQCD with full QCD, 
while yUf ac t denotes the normalization point for the calculations in the perturbation theory of 
NRQCD. 

For the heavy quarkonium composed of heavy quarks with the same flavor, the Wilson coeffi- 
cient /C is known with the two-loop accuracy fl3|-|T6l| 



«f (/ihard) , /af(/i hard r 2 



£■ (/ihard! /ifact) = 1 - ~~ ^ ~ c 2 (/ihard ', /ifact) , (2) 

O 7T \ 7T 



and C2 is explicitly given in [14],[15[]. The anomalous dimension of factor /C(/ih ar d; /ifact) in NRQCD 
is defined by 

dln/C(/z hard ;/i) -A f af s (fi)\ 

— — = ^ k] [-^r) ' (3) 

whereas the two-loop calculations^ give 

T[i] = 0, (4) 
7[2] = -16n 2 C F (^C F + ^C A y (5) 

The initial condition for the evolution of factor /C(/ihard! /ifact) is given by the matching of NRQCD 
current with full QCD at /i = /ihard [|14l,|T5j. 

The leptonic constant of heavy quarkonium is defined in the following way: 

(0\J? CD \QQ,\) = e x J QQ M QQl (6) 



x We use ordinary notations for the invariants of SU(N C ) representations: Cf — N 2 C N 1 , Ca — N c , Tf — 5, n/ is 
a number of "active" light quark flavors. 
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where A denotes the vector state polarization e v . In full QCD the electromagnetic current of quarks 
is conserved, while in NRQCD the current J^RQCD ^ e nonzero anomalous dimension, so that 
in accordance with (|1])-(|5|), we find 



(0\J? RQCD (v)\QQ, A) = AM e}f™ QCD M ( 



L QQ> 



(7) 



where, in terms of nonrelativistic quarks, the leptonic constant for the heavy quarkonium is given 
by the well-known relation with the wave function at the origin 



NRQCD 
QQ 



,12 , 



(8) 



and the value of wave function in the leading order is determined by the solution of Schrodinger 
equation with the static potential, so that we isolate the scale dependence of NRQCD current in 
the factor A(fi), while the leptonic constant Iqq® 00 is evaluated at a fixed normalization point 
fj, = Ho, which will be attributed below. It is evident that 



Jqq — /, 



NRQCD 



^(/^fact) • ^(/ihard; /^fact) 



(9) 



and the anomalous dimension of A(fifact) should compensate the anomalous dimension of factor 
^(/^hard; jUfect), so that in two loops we have got 



dliaA(fi) 
<iln/i 



-7[2] 



a: 



4tt 



(10) 



The physical meaning of A(/a) is clearly determined by the relations of (0) and (|j): this factor gives 
the normalization of matrix element for the current of nonrelativistic quarks expressed in terms of 
wave function for the two-particle quark state (in the leading order of inverse heavy quark mass in 
NRQCD). Certainly, in this approach the current of nonrelativistic quarks is factorized from the 
quark-gluon sea, which is a necessary attribute of hadronic state, so that, in general, this physical 
state can be only approximately represented as the two-quark bound state. In the consideration 
of leptonic constants in the framework of NRQCD, this approximation requires to introduce the 
normalization factor A(p) depending on the scale. 



The renormalization group equation of fllCf ) is simply integrated out, so that 



f3o + (3: 



47T 



A) + /3: 



oZ*<JH>) 

47T 



7[2] 

2fa 



(11) 



where (3q = —Ca — %TpUf, and f3\ = yCj — ACpTpUf — ^(TaTfti/. A constant of integration 
could be defined so that at a scale fio we would get A(no) — 1- Thus, in the framework of NRQCD 
we have got the parametric dependence of leptonic constant estimates on the scale fio, which 
has the following simple interpretation: the normalization of matrix element for the current of 
nonrelativistic quarks at /^o is completely given by the wave function of two-quark bound state. 
At other fi ^ fio we have to introduce the factor A(fi) ^ 1, so that the approximation of hadronic 
state by the two-quark one becomes inexact. 
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The relations above are general for the NRQCD description of leptonic constant for the heavy 
quarkonium in both the approach of sum rules and the calculations with the static potential, so 
that the only difference is the use of (§) in the potential approach. In the sum rules this relation is 
substituted by the equations for the constants I'qq® 00 , which follow from the quark-hadron duality 
in the corresponding scheme of calculations. In this respect, we have to stress that the analysis 
of NRQCD sum rules for the bottomonium in the two-loop approximation in || [7|, [16| took into 
account the factor of A(n) in a different way. So, in the calculation of two-point correlator for the 
currents of nonrelativistic quarks the regularization is necessary in the two-loop order. Of course, 
the anomalous dimension of two-point current correlator compensates the anomalous dimension of 
corresponding Wilson coefficient. In the sum rules with the one-loop gluon correction M this factor 
can be omitted, since its one-loop anomalous dimension is equal to zero, and, hence, the factor 
can be supposed equal to unit. As we have shown taking into account the two-loop corrections, 
the calculations without the normalization factor A(fi) for the matrix element of NRQCD current 
would be, in general, not correct. 

Some comments are to the point. It is worth to note that we deal with the operator product 
expansion (OPE), wherein the electromagnetic current of heavy quarks is expressed in terms of 
leading 1/mQ-order current of nonrelativistic spinors near the threshold. Therefore, in this OPE 
the corresponding Wilson coefficient, i.e. the matching factor /C, appears. The first problem is a 
calculation of this factor. The second is an evaluation of matrix element for the NRQCD current of 
nonrelativistic quarks. These two problems should be clearly distinguished in the estimates. The 
calculation of Wislon coefficient fC can be done perturbatively in two ways. Originally, it was done 
in the technique of threshold expansion of loop integrals for Feynman diagrams in full QCD for 
heavy quarks fl3|-|T6| . We use this result throughout this paper. Further, the second way of /C cal- 
culation is the following: one should perturbatively expand the QCD Lagrangian for the interacting 
heavy quark and heavy anti-quark separated by a distance r in terms of nonrelativistic spinors 
and represent the perturbative potential contributions different from the interactions with soft and 
ultrasoft fields in the hadron. These soft terms are not essential in the perturbative calculations of 
matching factor /C at high virtualities, and under such conditions the system of heavy quark and 
heavy anti-quark can be considered as the coulomb system in the leading approximation, while the 
potential corrections such as the spin- dependent and relativistic terms in the form of I/tuq con- 
tributions converted to the a s corrections to the energy, since in the coulomb systems the distance 
behaves as r ~ l/a s mg ~ l/|p| . Anyway, the main point is that the factor /C can be calculated 
in the perturbative theory. It is important that the fact of cancellation of anomalous dimension for 
the factor K, by the anomalous dimension for the NRQCD current follows from the zero anomalous 
dimension of electromagnetic current in full QCD with no reference to the calculation procedure, 
of course. Thus, the cancellation has no connection to the problem of evaluation for the hadronic 
matrix element of NRQCD current. Indeed, the effective expansion of QCD Lagrangian for the 
heavy quarks in the form of potential nonrelativistic QCD (pNRQCD) [[lTj or velocity-counting 
nonrelativistic QCD (vNRQCD, discussed below) |T8|-f20|1 cannot be applied to such the calculation 
of matrix element because the leading coulomb approximation does not dominate in comparison 
with corrections for the real quarkonia, even for the bb(lS). So, in pNRQCD and vNRQCD the 
potentials are treated as Wilson coefficients in front of four quark operators, while the soft and 
ultrasoft terms are separated. This separation is valid only if the time for the formation of coulomb 
wave function is many times less than the time of interaction with the soft fields (a quark-gluon 
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sea and a string). Then, 

T C oui T so ji. 

These times can be easily evaluated in terms of distance between the quarks, their velocity and 
soft scale Aqcd, so that 



r 

T cou l ~ - r ^ 



v rriQV 2 



rriQV 2 3> A 



QCD- 



T 



soft 



A 



QCD 



Therefore, in order to allow one to correctly use the pNRQCD or vNRQCD to the description of 
matrix elements, the kinetic energy of heavy quark motion should be many times greater than the 
scale of confinement. 

In practice, estimates for the most deep 'coulomb' system 66(15) give 



m Q v z ~ 500 MeV, at A QCD ~ 300 MeV. 

Thus, in order to get the potential model applicable to the real quarkonia, one should include soft 
terms in the static energy, i.e. the static potential. This static potential has an infrared stability 
and it is independent of any separation of soft and non-soft terms in the static energy. The static 
potential defined in terms of Wilson loop possesses these properties and it does not coincide with 
the potential (Wilson coefficient) in pNRQCD or vNRQCD, since the infrared contributions are 
added, but separated. Therefore, we deal with the model, wherein the total energy of all dynamical 
fields in the presence of static sources is taken into account. By the way, the relativistic and spin- 
dependent corrections in this case determine the scale of corrections about 50 MeV and their 
contribution gives the uncertainty of estimation, which is many times less than the uncertainty 
due to the scale variation in the estimates for the leptonic constants with the static potential. 

We stress that there is an ordinary situation in the leptonic constant calculation, since the 
perturbative technique used to evaluate the Wislon coefficient cannot be explored in order to get a 
correct estimate for the hadronic matrix element for the operator in front of the Wislon coefficient. 

Recently, the effective theory for the heavy nonrelativistic quark and antiquark pair in QCD 



was developed on the basis of relative velocity expansion [|T8|-|20"| . This approach uses the so-called 
velocity renormalization group, so that we refer to such the theory as vNRQCD. The physical 
argumentation of vNRQCD is the following: In the heavy quark-antiquark system one expects a 
low influence of light quark-gluon sea formed by nonperturbative emission of quarks and gluons with 
virtualities about Aqcd, so that the heavy quark motion is close to the coulomb approximation up 
to higher order corrections in a s . Moreover, we can suppose the coulomb relations between various 
scales, which possess the following counting rules: the relative velocity v ~ a s , the momentum 
transfer (or the inverse distance between the quarks), i.e. the soft scale fis ~ |fe| ~ mg • v, the 
momentum of heavy quark |p| ~ fi$, while the kinetic and potential energies, i.e. the ultrasoft 
scale, [Au ~ E ~ V ~ uiq ■ v 2 , so that in the perturbative part of effective Lagrangian the distance, 
momentum and energy are not free independent quantities, since they are related by the single 
"running" parameter v. This approach is certainly different from pNRQCD, wherein the static 
bilinear operators of nonrelativistic spinors are fixed at a short distance, while the kinetic energy 
as well as the quark momenta are put to zero. 

Thus, the starting point of vNRQCD is the correlation of soft and ultrasoft scales, nu « v ■ us, 
in contrast to pNRQCD. Therefore, in vNRQCD the matching with the full QCD at the scale of 
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heavy quark mass /Xhard — and the renormalization group evolution with v involving large In v 
result in a combined summationP] of relevant momentum and energy contributions from fi = itlq 
up to n — rriQ ■ v 2 , while pNRQCD deals with the two-step matching and evolution: NRQCD 
matched with full QCD at /ihard ~ m Q and evolved to /is ~ rriQ ■ v, and pNRQCD matched with 
NRQCD at fi s and evolved to fi ~ mg ■ v 2 . A peculiar feature of vNRQCD summing up both 
logs of fis and \iu is that the limits of mq — > oo or v — > do not reproduce the static limit (see 
explanations provided by examples in |20(|). 



Next point is the anomalous dimension of nonrelativistic current, that was calculated in [p!q ,p!9|] 
in the technique of vNRQCD. One found that expressions are reproduced in vNRQCD 

and strictly valid at /i = tuq, while at lower virtualities with v = ^f ac t/ m Q < 1 some corrections 
in \nv appear. These corrections can be summed up in vNRQCD, so that substitution of a 2 {^) 
for a 2 [mo) in the two- loop anomalous dimension could be quite naive and not strictly justified, 
by opinion in ||19|| . The analogous two- loop anomalous dimension 7 of K, in vNRQCD as well as 



the complete expression for /C vNR Q° D (t>) integrated out in the velocity renormalization group are 
available in [Ofl. So, 



^vNRQCD (|;) = £vNRQCD (1) exp 



S ^2a v i v:a s (mQ)ki 



i=l 



(12) 



where are expressed in terms of Cp, Ca and /3q, and the functions ki have the form 

k\ = - - 1, k-2 = 1 — z, k 3 = ln(A 

z 

jfe 4 = 1 - Z l-13C A /(6/3 ) ) jfe 5 = 1 _ 2 l-2C A //3o 5 

k( , = t. _ I i n 2 (2 ) _ l n ( w ) l n f ^^^) _ Li 2 ( — 



12 2 v ' \2w -I J \2w 

w 1 
k 7 = 2w - 1 ln(iw) - - ln(2tu - 1), 

while z = a s (/i)/« s (mg), w = a s (/x 2 / mg) / a s (/x) , and the one-loop matching condition is given by 
the same form as in (2), i.e. 



a 



: (m Q ) 



^viNuycu/^ = i _2C F (13) 

7T 

However, in general the scale prescription in the argument of a s standing in the first nonzero term, 
which presents the accuracy under study, is rather conventional, since the scale change results in 
the next-order contribution in a s , that is not under control. Indeed, the substitution 



a s (yu) ~ a s (m Q ) ( 1 - ^ a s (m Q ) ln — ) 
V. 2vr m Q J 



a s (m Q ) ( 1 - ^ a s (m Q ) Inv 
In 



influences the three-loop anomalous dimension 7 for /C, only. The argumentation of vNRQCD is 
that the corrections summing up [lnf] n are dominant (in all orders in a s ). So, the above note on 
the scale substitution manifests itself in three-loop running of coulomb-like potentials in vNRQCD. 

In the present paper we follow the static ideology with the threshold expression for fC(fi) (and 
hence, A(fi)), since we adopt the Schrodinger equation with the model of static potential for the 

2 Note that terms [lnu] ra ~ [lna s ] ra reproduce so-called manifestly nonanalytic contributions in o^-expansion. 
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QQ system in contrast to the coulomb approximation improved by the velocity renormalization 
group involving running coulomb and relevant potentials. Moreover, the threshold expression for 
/C(/a) is preferable by technical reasons, since, first, it is valid at the two- loop order (as we have 
mentioned above), second, /C(/i) has an explicit dependence on the matching scale /ihard, that 
allows us to study the point stable under the variation of /ihard (see below), while in vNRQCD 
the matching scale is fixed and prescribed to /ihard = m Q, and the corresponding analysis on the 
stability is not possible. Third, the factor /C(/i) contains the two-loop final term independent of 
In /i, while /C vNR< ^ CD (v) does not (except the one-loop term), since in the large lnw approximation 
the two-loop final term is irrelevant. This two-loop final term in /C(/i) makes the cancellation of 
/i-dependence in the product of /C(/i) • A(n) valid in a restricted region of /i, that is conceptually 
true, while in vNRQCD the large log approximation results in the complete cancellation of in- 
dependence (remember, v = /ifact/^g) in the product of /C vNR Q CD (/i) ■ ^4 vNR< ^ CD (/x). Nevertheless, 
since we assume the existence of /io, at which A(fio) = 1, then we can estimate the above product 
by /C vNR( ^ CD (wo). The problem on the choice of /io in vNRQCD is more artificial than in NRQCD, 
since the analysis of stability on /ihard results in a preferable prescription for /i in NRQCD, while 
in vNRQCD we have no analogous criterion (see Section 4 with numerical results). 



3 Potential approach 



The potential of static heavy quarks illuminates the most important features of QCD dynamics: the 
asymptotic freedom and confinement. In the leading order of perturbative QCD at short distances 
and with a linear confining term in the infrared region, the potential of static heavy quarks was 
considered in the Cornell model 
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incorporating the simple superposition of both asymptotic 
limits (the effective coulomb and string-like interactions). The observed heavy quarkonia posed in 
the intermediate distances, where both terms are important for the determination of mass spectra 
(see Fig. [I]). So, the phenomenological approximations of potential (logarithmic one [^] and power 
taking into account the regularities of such the spectra, were quite successful p|[ . 



law |f23 



The quantities more sensitive to the global properties of potential are the wave functions at 
the origin as related to the leptonic constants and production rates. So, the potentials consistent 
with the asymptotic freedom to one and two loops as well as the linear confinement were proposed 
by Richardson f25|, Buchmiiller and Tye |L2|, respectively. 



In QCD the static potential is defined in a manifestly gauge invariant way by means of the 
vacuum expectation value of a Wilson loop [p6fl, 



V(r) 



Wr 



-Km -ln(Wr) 

T^oo ll 



tr V exp 



(14) 



Here, T is taken as a rectangular loop with time extension T and spatial extension r. The gauge 
fields are path-ordered along the loop, while the color trace is normalized according to tr(..) = 
tr(..)/trl. This definition corresponds to the calculation of effective action for the case of two 
external sources fixed at a distance r during an infinitely long time period T, so that the time- 
ordering coincides with the path-ordering. Moreover, the contribution into the effective action by 
the path parts, where the charges have been separated to the finite distance during a finite time, 
can be neglected in comparison with the infinitely growing term of V(r) - T. Let us emphasize that 
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Figure 1: The Cornell model of static potential and sizes of observed heavy quarkonia with charmed 
quarks (the family ip) and bottom quarks (the family T). 

the defined static potential is, by construction, the renormalization invariant quantity, since the 
action, by definition, doe not depend on the normalization point. 

Generally, one introduces the V scheme of QCD coupling constant by the definition of QCD 
potential of static quarks in momentum space as follows: 

= _c/Jl^l, (15) 

so that for the such-way introduced value of a v one can derive some results at large virtualities in 
the perturbative QCD as well as at low transfer momenta in the approximation of linear term in 
the potential confining the quarks. 

In this section, first, we discuss two regimes for the QCD forces between the static heavy quarks: 
the asymptotic freedom and confinement. Then we follow the method by Buchmiiller and Tye and 
formulate how these regimes can be combined in a unified f3 function for a v obeyed both limits of 
small and large QCD couplings. 

In the perturbative QCD up to the two-loop accuracy, the quantity a v can be matched with 
«ms by the relation 
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At present, in (|T6| ) the coefficients of tree approximation a , one-loop contribution a% and new 
two-loop results a 2 (see |IO|,|ll| ) are known. Note that expansion (|16D cannot be straightforwardly 
extended to higher orders of perturbative QCD because of infrared problems, that result in non- 
analytic terms in the three-loop perturbative potential as was first discussed by Appelquist, Dine 
and Muzinich 



After the introduction of a = j-, the (3 function is actually defined by 

da(fi 2 



din n 2 



/3(a) = -]T/3„-a" + V), (17) 



n=0 



so that f% fl = 01$ and /3 2 v = (3f s - a x (3f B + (a 2 - aj)f3™ s . 

The Fourier transform results in the position-space potential [10|, so that the perturbative 
potential, by construction, is independent of normalization point, i.e. it is the renormalization 
group invariant. However, in the problem under consideration the truncation of perturbative 
expansion, wherein the coefficients do not decrease^, leads to a strong custodial dependence on 
the normalization point. So, putting the normalization point \i in the region of charmed quark 
mass, we find that the two-loop potential with the three-loop running coupling constant af s has 
an unremovable additive shift depending on fi. This shift has variation in wide limits. This 
fact illuminates the presence of infrared singularity in the coupling constant of QCD, so that the 
//-dependent shift in the potential energy has the form of a pole posed at Aqcd 0- 

Thus, in order to avoid the ambiguity of static potential in QCD we have to deal with infrared 
stable quantities. The motivation by Buchmuller and Tye was to write down the (3 function of 
a v consistent with two known asymptotic regimes at short and long distances. They proposed 
the function, which results in the effective charge determined by two parameters, only: the per- 
turbative parameter is the scale in the running of coupling constant at large virtualities and the 
nonperturbative parameter is the string tension. The necessary inputs are the coefficients of (3 
function. The parameters of potential by Buchmuller and Tye were fixed by fitting the mass spec- 
tra of charmonium and bottomonium 0. Particularly, in such the phenomenological approach 



the scale A"p « 510 MeV was determined. It determines the asymptotic behaviour of coupling 
constant at large virtualities in QCD. This value is in a deep contradiction with the current data on 
the QCD coupling constant af s ||. In addition, one can easily find that the three-loop coefficient 
P2 for the (3 function suggested by Buchmuller and Tye is not correct even by its sign and absolute 
value in comparison with the exact coefficient recently calculated in ]T0|,|n{]. 

Thus, the modification of Buchmuller-Tye (BT) potential of static quarks as dictated by the 
current status of perturbative calculations is of great interest. 

The nonperturbative behaviour of QCD forces between the static heavy quarks at long distances 
r is usually represented by the linear potential (see discussion in ref. |27| ) 

V conl {r) = k-r, (18) 

which corresponds to the square-law limit for the Wilson loop. The form of ( |T8"D corresponds to 
the limit, when at low virtualities q 2 — > the coupling a v tends to 

"v(q ) -> -j, 

q 



3 Moreover, according to the investigations of renormalon, the coefficients in the series of perturbation theory for 
the potential increase in the factorial power, so that the series has a meaning of asymptotic one. 
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so that 



da v (q 2 



-a v {q 



(19) 



din q 2 

which gives the confinement asymptotics for the /3 V function. 

Buchmiiller and Tye proposed the procedure for the reconstruction of (3 function in the whole 
region of charge variation by the known limits of asymptotic freedom to a given order in a s 
and confinement regime. Generalizing their method, the /?pt function found in the framework of 
asymptotic perturbative theory (PT) to three loops, is transformed to the (3 function of effective 
charge as follows: 



1 



Apt (a) 

1 



(3(a) 



P a 2 ( 1 — exp 



i 



+ 



A + (g ~ |) 
P 2 a 



exp 



l 2 a 2 



(20) 



where the exponential factor in the second term contributes to the next-to-next-to-leading order 
at a — > 0. This function has the essential peculiarity at a — > 0, so that the expansion is the 
asymptotic series in a. At a — > oo the f3 function tends to the confinement limit represented in 

The construction of (p0[) is based on the idea to remove the pole from the coupling constant 
at a finite energy, but in contrast to the "analytic" approach developed in [^] , the "smothing" of 
peculiarity occurs in the logarithmic derivative of charge related with the f3 function, but in the 
expression for the charge itself. Indeed, in the one-loop perturbation theory we have got 



din a 
dln/i 2 



because of 



-A) a 



1 



ln^' 



and the pole can be cancelled in the logarithmic derivative itself, so that[] 



din a 

din /i 2 



A 2 

-fta(l- ? 



(3oa ( 1 — exp 



As we see in the perturbative limit o — > the deviation in the f3 function is exponentially small. 
The generalization of this idea to two loops was done by Buchmiiller and Tye. 

Remember that the one and two-loop static potentials matched with the linear term of confine- 
ment lead to the contradiction with the value of QCD coupling constant extracted at the scale of 
Z boson mass if we fit the mass spectra of heavy quarkonia in such potentials. The static potential 
in the three-loop approximation results in the consistent value of QCD coupling constant at large 
virtualities M. 



'The presentation of power correction in the form of charge exponential was similarly used in p£ 
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In the perturbative limit the usual solution for the running coupling constant 



a(/i 5 



1 



i & 1 i i »\ 
1 - — in In — + 



A 2 



01 



/?oln 2 ^ 



In In — — In In — — 1 + 



A 2 



A 2 



A 2 



is valid. Using the asymptotic limit of (|21f) , one can get the equation 

u 2 1 /3i /" a(At2) 

A 2 A)a(/i 2 ) /?o 7o 



A + 1 



fox 2 (3 2 x (3(x) 



(21) 



(22) 



which can be easily integrated out, so that we get an implicit solution for the charge depending 
on the scale. The implicit equation can be inverted by the iteration procedure, so that well 
approximated solution has the form 



a(/i 2 



A)ln(l+i?(A* a )fr 



(23) 



where r](fi 2 ) is expressed through the coefficients of perturbative (3 function and parameter I in 
(|20|), which is related to the slope of Regge trajectories a' P and the integration constant, the scale 
A@. 

The slope of Regge trajectories, determining the linear part of potential, is supposed equal to 



1.04 GeV , so that in (18) we put the parameter k 



2na' 



We use also the measured value 



of QCD coupling constant M and pose 



offal) = 0.123, 

as the basic input of the potential. The transformation into the configuration space was done 
numerically in 0, so that the potential is presented in the form of file in the notebook format of 
MATHEMATICA system. 

The analysis of quark masses and mass spectra of heavy quarkonia results in the following 
values ascribed to the potential approach M: 



1.468 GeV, ml 



4.873 GeV. 



(24) 



Thus, the spectroscopic characteristics of systems composed of nonrelativistic heavy quarks are 
determined in the approach with the static potential described above. 



4 Numerical estimates 

In the potential approach the values of wave functions for the IS" states in the systems of bb and 
cc are equal to 

>/&r a (0) = 2.513 GeV 3/2 , \/4^ 5c (0) = 0.895 GeV 3/2 , (25) 

which are the renormalization group invariants by construction of Schrodinger equation with the 
static potential. Remember that this potential includes the contribution of linear part providing 
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the confinement of quarks. The same term makes the infrared stability of effective charge and all 
of the potential in NRQCD. 

The analysis of leptonic constants of T and ip with account of Wilson coefficient /C, derived in 
the perturbative NRQCD, was done in ||. For T we found that the variation of hard scale in wide 
limits /ihard = (1 — 2)mb led to the presence of a stable point /if act , in which the result was weakly 
sensitive to the variation of hard scale. The stability was observed at /i fact pa 2.2 — 2.7 GeV, i.e. 
at the inter-quark distances characteristic for the size of IS level in the system of bb. 

For the leptonic constant of charmonium J/ip the stability point under the variation of /ih ar d 
was located at the reasonable values of soft scale of factorization /if act pa 1.05 — 1.35 GeV. The 
position of stable point was significantly dependent of changes in /ihard close to the mass of charmed 
quark. 

However, in that consideration the stability of results under the variation of scale /if act for the 
perturbative calculations concerning the Wislon coefficient in NRQCD was not observed. 

In this paper we introduce the normalization factor A into the analysis for the current of 
nonrelativistic quarks. As we see in Fig. |2], in the system of bb the choice of initial point of 
normalization /a = 2.3 GeV in the mentioned region of stability for the leptonic constant under 
the variation of matching scale, so that .4.(2.3 GeV) = 1, we get the independence of result for the 
lepton constant of T in the interval /if act = 1.8 — 3.8 GeV. 

fr, GeV 

0.8 



0.7 



0.6 



0.5 



!-5 2 2.5 3 3.5 4 GeV 



Figure 2: The leptonic constant of ground vector state in the system of bottomonium is presented 
versus the soft scale of normalization. The dashed curve corresponds to /ihard = 2mb, while the 
solid curve does to /ih ar d — m b- The initial condition for the evolution of normalization factor 
A(fif act ) in the matrix element of current in the nonrelativistic representation has been posed in 
the form ^4(2.3 GeV) = 1. The horizontal band is the experimental limits for the constant. 

The analysis of stability is performed in Fig. ^. Thus, our estimate of leptonic constant for the 
ground vector state of bottomonium is given by 

fr = 685 ± 30 MeV, 
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/i fact , GeV 



/ihard, GeV 



Figure 3: The leptonic constant of ground vector state in the system of bottomonium is presented 
versus the soft and hard scales, /if act and /ihard- The initial condition for the evolution of nor- 
malization factor 4.(/if ac t) in the matrix element of current in the nonrelativistic representation is 
.4(2.3 GeV) = 1. 



GeV 

0.6 



0.5 




/i fact , GeV 



Figure 4: The leptonic constant of ground vector state in the system of charmonium is presented 
versus the soft scale of normalization. The shaded region restricted by curves corresponds to the 
change of hard scale from /ihard = 2.0 m c (the dashed curve) to /ihard = 1-3 m c (the solid curve). 
The initial condition for the evolution of normalization factor 4(/if ac t) in the matrix element of 
current in the nonrelativistic representation has been posed in the form 4.(1.07 GeV) = 1. The 
horizontal band is the experimental limits for the constant. The additional curves are presented 
for /ihard — 2.8 m c (the dashed line) and /ihard = 1.15 m c (the solid line). 
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which is in a good agreement with the experimental data. 

The calculations of leptonic constant for the charmonium J / ip are more sensitive to the choice 
of initial condition in the evolution of factor A, because the stability point under the variation 
of hard scale for the matching is significantly fluctuates under the change of /ihard = (1 — 2) m c . 
In Fig. [| we show the dependence of fu^ on the scale //f act at .4,(1.07 GeV) = 1. This value of 
fiQ = 1.07 GeV agrees with the inverse size of system cc. We see in the figure that the region of 
result stability is given by //f act = 1. - 1.6 GeV. 

/j/^, GeV 




Figure 5: The leptonic constant of ground vector state in the system of charmonium is presented 
versus the soft and hard scales, /if act and ^hard- The initial condition for the evolution of normaliza- 
tion factor *4.(/ifact) in the matrix element of current in the nonrelativistic representation is given 
by .4(1.07 GeV) = 1. 

The dependence of estimates on the scales are presented in Fig. [5[ Then we get 

fj/4, = 400 ± 45 MeV, 

which can be compared with the experimental value /J^ = 409 ± 15 MeV. 

Thus, the analysis of stability for the leptonic constants of heavy quarkonia in the framework 
of potential approach shows that there is the region of stability under the change of both factor- 
ization scale in the perturbative NRQCD and point marking the matching of effective theory of 
nonrelativistic quarks with the full QCD. The basic source of uncertainty in the estimates is the 
variation of initial normalization point for the evolution factor determining the matrix element for 
the current of nonrelativistic quarks, so that the values of leptonic constants in the extremum (see 
Figs. ^| and |5[) change in the limits, which are included into the methodic error of estimates as 
shown above. 

In vNRQCD the factor /C vNRQCE » demonstrates the stability region (see Fig. || for the bot- 
tomonium). This stable point is the only reasonable indication for the choice of vq, at which 
*4 vNR Q CD (v ) = 1. Such the prescription results in the estimates of leptonic constants derived in 
vNRQCD 

/r| vNRQCD = 680 ±30MeV, f m | D = 360 ± 45 MeV, (26) 
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where the errors are dominantly due to the variations of a s at the masses of bottom and charmed 
quarks: a s (mb) = 0.22 ±0.02, a s {m c ) = 0.34 ±0.05 entering the one-loop matching condition. We 
see that the vNRQCD estimates are in agreement with both the static values and experimental 
data. 

£vNRQCD^ 
£vNRQCD(]_) 

1.1 
1.05 
1 

0.95 
0.9 

2.5 3 3.5 4 4.5 

ji = niQ ■ v, GeV 

Figure 6: The Wilson coefficient /C vNRC ^ CD (v) for the vector electromagnetic current of nonrela- 
tivistic bottom quarks in vNRQCD. 

Finally, in the potential approach we can calculate the ratios of leptonic constants for the excited 
nS-waye levels in bb and cc, which were presented in U in comparison with the experimental data. 
These theoretical calculations are in a good agreement with the measured values. 

5 Conclusion 

In the present work we have done the general consideration of problem on the calculation of leptonic 
constants for the heavy quarkonia in NRQCD and introduced the factor for the normalization 
of matrix element for the current of nonrelativistic quarks. This factor has got the anomalous 
dimension of current in the effective theory, so that the explicit invariance under the action of 
renormalization group takes place for the estimated values of leptonic constants. In the framework 
of potential approach we have estimated the constants for /x and fj/$, so that we have found 
reasonable regions of stability under the variation of both the scale for the matching of NRQCD 
with full QCD and the soft scale ascribed to the operators of nonrelativistic quarks in NRQCD. 
The obtained results are in a good agreement with the experimental data. 
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